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1.  INTRODUCTfON 


nio  now  of  a  liquid  jet  inlo  a  near  vacuum  envirenment  has  lieen  studied  by  many  investigators.  Many 
of  llifsc  studies  are  experimental  but  also  provide  some  theoretical  interpretation  of  the  experiments, 
l  or  exam()le,  see  the  works  of  Joltnson  ;md  Hullett  |l'X)Kj.1C'issal  1 1974|,  Fuchs  and  Legge  |1979],  Curry 
et  al.  II9S5a,b|,  and  Muntz  and  Orme  II987J.  Johnson  and  Ilallelt  |1%8]  investigated  the  behavior  of 
supera)oled  water  drops  at  various  conditions,  such  as  in  carbon  dioxide  below  0.05  atmospheric 
pressure.  Kassal’s  ( 1974)  experiment  showed  that  the  spherical  particles  break  in  half  during  the  freezing 
process  and  form  bowl-shaped  configurations  in  the  vacuum  chamber.  Fuchs  and  Legge  [1979] 
concluded  that  the  bursting  effect  caused  by  the  sudden  boiling  of  the  water  produces  vapor  bubbles 
which  decompose  the  jet.  They  stated  that  the  formation  of  vapor  bubbles  was  caused  by  a  suf>erheating 
of  the  water  due  to  sudden  pressure  drop  and  cooling  at  the  surface  of  the  jet.  Muntz  and  Orme  [1987] 
discussed  the  possible  events  when  liquid  streams  are  released  into  a  vacuum.  'Iliey  stated  that  for  fluids 
having  vapor  pressures  up  to  several  torr  in  streams  with  a  few  hundred  microns  in  diameter,  it  is 
possible  to  maintain  stable  and  controllable  streams.  However,  fluid  streams  having  vapor  pressures  in 
the  of  tens  of  torr  range  with  one  millimeter  and  greater  diameter  are  expected  to  burst. 

Curry  et  al.  [1985a,  b]  performed  an  experiment  in  a  vacuum  chamber  (10'®  Torr),  wherein  water  was 
released  from  an  orifice  with  a  diameter  of  1.6  mm  into  the  chamber.  TTie  vent  plume  was  studied  using 
a  Mie  scattering  polarimeter  and  a  high-speed  flow  visualization  camera.  They  concluded  that  there  are 
mainly  two  classes  of  particles  which  can  be  found  in  the  water  plume,  viz.  the  large  particles  with 
diameters  close  to  the  orifice  diameter  (1.6  mm)  and  a  cloud  of  very  small  particles  whose  diameters  were 
about  0.15  pm.  Pike  et  al.  [1990]  showed  that  the  particle  sizes  as  measured  from  the  recent 
groundbased  and  onboard  video  images  of  a  sunlit  Shuttle  Orbiter  water  dump  are  similar  to  those  noted 
in  vacuum-tank  experiments  (Curry  et  al.  [1985a,b]). 

Some  mathematical  modeling  of  such  flows  has  also  been  done  and  compared  with  experimental  data 
(Gale  et  al.,  1964,  Glenn  1969).  Gale  et  al.  [1964],  in  their  mathematical  study  of  liquids  (including 
water)  released  into  a  vacuum,  analyzed  the  effects  of  initial  temperature,  droplet  diameter,  and  ambient 
pressure  on  the  time  required  for  cooling  to  the  triple  point  and  for  subsequent  freezing.  Part  of  their 
mathematical  analysis  was  supported  by  experimental  data.  Glenn  [1969]  studied  a  radial  source  flow 
of  uniform-size  liquid  particles  into  a  vacuum.  A  numerical  procedure  for  the  integration  was  developed 
for  the  models,  and  he  discussed  parametrically  the  effects  on  the  flow  structure  of  such  variables  as 
particle  size,  velocity,  and  temperature,  llie  present  investigation  is  inspired  by  the  work  of  Glenn  [1969] 
in  which  h(  assumed  a  one-dimensional  jet  flow  of  spherical  particles  of  liquid  water  into  vacuum. 
Glenn’s  [1%9]  work  has  been  extended,  and  some  revisions  made  to  his  work  which  are  included  in  this 
report. 


2.  DESCRIPTION  OF  THE  MODEL 


I  he  problem  considered  here  is  the  mathematical  modeling  of  stationaiy  radial  source  flow  into  a 
vacuum.  Ibe  equations  describing  this  model  are  the  one-dimensional  continuity,  momentum,  and 
energy  for  the  vapor  and  particles.  The  analysis  constitutes  a  sixth-order,  nonlinear,  two-point  boundary 
value  problem.  A  numl>er  of  assumptions  have  been  made  to  simplify  the  problem.  These  include  the 
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following:  'llie  particles  are  liquid,  spherical  shaped,  with  uniform  size  and  temperature.  Hie  wake  of 
the  particles  does  not  influence  the  nearby  panicles.  Re-condensation  of  vapor  as  a  result  of  any 
subcooling  during  the  expansion  process  is  neglected.  ITie  heat  transfer  between  particles  and  vapor 
occurs  only  by  convection  and  radiation.  Ihe  fluid  is  considered  to  be  inviscid,  except  the  part  that  can 
produce  drag.  From  the  above  constraints,  one  can  assume  a  simple  liquid  jet  breakup  into  the  vacuum 
environment. 

ITie  derivations  of  the  equations  that  describe  this  model  am  be  found  in  Appendix  A.  Glenn’s  work 
[1969]  has  been  extended  and  the  equations  have  been  modified  to  investigate  the  behavior  of  particles 
after  freezing,  and  to  include  heat  radiation  from  the  particle.  The  Oausius-Clapeyron  equation  is  used 
for  calculating  the  local  vapor  pressure  in  all  cases  considered.  The  problem  has  been  solved  with  the 
engineering  parameters  for  the  release  of  water  into  space.  ITiese  parameters  were  chosen  to  resemble 
those  ot  a  recent  water  release  experiment  from  the  space  shuttle  (Discovery,  March  1989).  Nevertheless, 
it  should  be  mentioned  here  that  the  mechanism  for  water  release  and  its  breakup  into  space  is  an 
extremely  complex  process  which  cannot  be  fully  described  by  this  simple  one-dimensional  model. 
However,  it  may  give  some  indication  of  the  factors  involved.  Future  work  will  include  two,  and  perhaps, 
three  dimensional  models. 

'Hie  governing  equations  that  are  used  in  this  study  are  as  follows: 
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I:qu;iti(>ns  ( 1  -6)  coiistituto  :i  sixth  order  system  in  six  unknown  quantities;  vapor  density  p^,  vapor  velocity 
11^.  particle  velocity  Up,  particle  temperature  I'p,  vapor  temperature  1’^,  and  particle  volume  fraction  «, 
and  five  additional  variables  which  will  be  defined  later. 

Ibe  above  equations  are  subject  to  the  following  constraints: 


r  >  T.  ;  p  =  1.  =  0 

T,  =  r.. ;  t;:  =  0 

Ty<T,,-,  p  =  p'  =  0 


Where  f..  and  'f.  are  the  triple  point  and  initial  temperatures  respectively,  and  p  represents  the  particle 
liquid  volume  fraction. 

Ilie  five  additional  unknowns,  mass  rate  m,  particle  radius  a,  local  pressure  P,  drag  force  fp,  and  total 
heat  tr.insfer  by  radiation  and  convection  qp,  are  defined  as  follows: 

m  =  A  [Po  (^p)  -  P]  (2tt  R  Ty)-  (8) 
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llie  following  qinintities,  based  on  the  initial  conditions,  can  be  introduced  in  the  problem: 
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nien  the  results  cm  be  written  in  matrix  form: 
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■Jlic  O  in  the  sixth  row  and  fourth  column  of  the  matrix  is  equal  to  0g,  except  at  the  triple  point  Tp  =» 
1'..,  where  0  =  0|q;  and  parameters  A,  thron»h  are  defined  as  follows: 
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Here,  the  loail  vapor  pressure  as  a  function  of  particle  temperature  can  be  stated  in  terms  of  the 
Q:uisius-Qapeyron  equation 


fW  =  * 


1 


(22) 


3.  METHOD  OF  SOLUTION 


The  finite  difference  and  shooting  methods  are  two  of  a  number  of  methods  available  in  the  literature 
to  solve  the  two-point  boundary  value  problems.  Here  the  shooting  method  was  chosen.  Glenn  [1968] 
used  the  simplest  form  of  shooting  method,  called  the  "  straight-shooting”  scheme.  In  this  method 
(which  is  also  common  in  all  other  shooting  methods),  the  two-point  boundary  value  problem  can  be 
reduced  to  an  initial  value  problem.  The  differential  equations  are  integrated  as  an  initial  value  problem, 
and  the  initial  conditions  are  adjusted  until  the  conditions  at  the  other  end  point  of  the  interval  are 
satisfied.  Further  discussion  of  this  method  can  be  found  in  Glenn’s  work  [1S)68]. 


4.  DISCUSSION  OF  THE  RESULTS 


Figures  1  through  5  present  comparisons  of  the  current  work  to  the  work  of  Glenn  [1969],  and  Fuchs 
and  Legge  [1979].  Figures  1  to  4  show  the  variation  of  the  particle  temperature,  particle  velocity,  vapor 
temperature,  and  particle  size  (evaporation  rate)  with  respect  to  radial  distance  r  from  the  source  for 
different  sizes  of  the  particles.  Most  of  the  initial  conditions  in  these  figures  are  similar  to  those  which 
were  used  by  Glenn  [1969].  For  example,  the  initial  temperature,  the  initial  velocity  of  the  stream,  and 
the  source  diameter  are  taken  to  be  93.65  ®C,  15.24  m/s,  and  6.35  mm,  respectively.  The  parameters  used 
for  the  purpose  of  comparison  with  the  work  of  Glenn  [1969],  and  Fuchs  and  Legge  [1979],  can  be  found 
in  Tables  1  and  2,  respectively.  Figure  1  compares  the  temperature  variation  for  different  initial  size 
particles.  In  Figure  2,  the  velocity  changes  are  investigated  and  compared  to  the  work  of  Glenn  [1969]. 
Since  the  vapor  velocity  increased  very  rapidly  after  injection  of  the  water  into  the  vacuum,  the  smaller 
particles  are  most  likely  affected  by  this  increase  in  vapor  speed  and  their  velocity  will  increase  according 
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Table  1.  Parameters  Used  for  Figures  1- 


Initial  I eni|K‘raiiiri.'  (°K) 

Initial  Velocity  (m/s) 

Initial  Source  Radius  (m) 

Initial  Particle  Radius  (jim) 

Initial  Vapor  Pressure  (N/m^) 

Initial  Vapor  Density  (kg/m^) 

Initial  V'ciume  Fraction 

Evaporation  Coefficient 

Drag  Coefficient 

Particle  Nusselt  Number 

V'apor  lliermal  Conductivity  (W/m  °K) 

Vapor  Specific  Heat  (J/kg  °K) 

Heat  of  Evaporation  (J/kg) 

Particle  Density  (kg/m^) 

Gas  Constant  for  Water  (J/Tcg  °K) 


3W).f)5 

15.24 

3.175  X  10'^ 

12.5,  25,  or  5.7 

8.0  X  1" 

0.47916 

0.74 

0.05 

0.9 

2.73 

0.024403 
1863.172 
2.2741  X  10® 
962.8 
461.94 
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Table  2.  Parameters  Used  for  Figure  5 


Initial  Temperature  (°K) 

Initial  Velocity  (m/s) 

Initial  Source  Radius  (m) 

Initial  Particle  Radius  (pm) 

Initial  Vapor  Pressure  (N/m^) 

Initial  Vapor  Density  (kg/m^) 

Initial  Volume  Fraction 

Evaporation  Coefficient 

Drag  Coefficient 

Particle  Nusselt  Number 

Vapor  Thermal  Condu( ..  .ty  (W/tn  °K) 

Vapor  Specific  Heat  (J/kg  "K) 

Heat  of  Evaporation  (J/kg) 

Heal  of  Fusion  (J/kg) 

Particle  Density  (kg'nt^) 

Gas  Constant  for  Water  (J/kg  ®K) 


300.15 

12 

1.0  X  lO"* 

5 

4.0  X  10^ 

U.02S 

0.74 

0.05 

i.j 

2.73 

0.024403 
1863.172 
2.386136  X  10* 
3.340777  X  10* 
987.1 


461.94 


ti»  Ilk*  p.irlicle  size.  For  cx:inii)lL‘,  (he  hirj^ier  prirlicles  will  he  less  iiffecleil  by  the  v;ipor  velocity.  'Iliere 
is  ;i  luriher  tliscussiuii  with  respect  to  the  velocity  increase  of  the  particles  in  (he  later  part  of  this  work. 
Figure  shows  (he  vapor  temperature  variation.  'Ilie  temperature  drop  in  vapor  is  .somewhat  different 
from  that  found  in  the  particles.  Figure  .1  shows  that  the  vapor  temperature  decreases  more  rapidly  for 
the  larger  particles  than  for  the  smaller  particle.s.  Figure  4  shows  the  size  reduction  of  the  particles  due 
to  evap(;ration,  as  a  function  of  the  radial  distance.  'ITii-s  figure  indicates  that  the  smaller  particles 
evaporate  much  faster  than  larger  particles  near  the  source.  As  am  be  seen  from  Figures  1  through  4, 
and  in  particular  Figure  2,  the  present  results  are  in  good  agreement  with  the  work  of  Glenn  |1969]. 

Figure  5  investigates  the  changes  in  vapor  pressure  as  a  function  of  radial  distance  R,  based  on  Equation 
(23).  'Ibis  equation  relates  the  radial  distance  R  to  the  Mach  number  and  the  specific  heat  ratio  of 
vapor,  rhe  r^^  is  the  fictitious  sonic  radius  defined  by  Fuchs  and  Legge  [1979].  Since  the  vapor  flow  is 
supersonic,  the  vapor  pressure  is  calculated  based  on  the  Rayleigh-Pitot  formula  (Equation  (24),  where 
Pg  is  defined  as  a  critical  pressure  at  sonic  r^^ )  which  is  considered  by  Fuchs  and  Legge  [1979J.  In  this 
figure  the  e.xperimental  measurement  of  Fuchs  and  Legge  [1979]  compares  ''ery  well  with  the  calculated 
results  of  the  present  work.  As  can  be  seen  from  Figure  5,  the  curves  start  at  some  distance  from  the 
source,  the  distance  th.it  the  vapor  needs  to  travel  to  reach  supersonic  speed. 
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2*{y-\)M^  2(Y-l) 
•Y+1 


(23) 
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Figures  6  through  12  show  water  flow  characteristics  near  the  source  of  release,  into  a  vacuum 
environment,  for  the  parameters  listed  in  Table  3. 

Figure  6  shows  the  variation  of  particle  temperature  as  a  function  of  initial  particle  size  for  three  different 
particles.  Ibe  particle  sizes  are  well  within  the  estimates  for  water  release  from  the  space  shuttle 
Discovery  for  the  flight  of  March  1989  (ranges  within  0.25  pm  to  1  millimeter  sizes).  The  initial 
temperature  of  the  liquid  water  is  70  "C  (.343.15  ®K),  which  leaves  with  an  initial  velocity  of  15  m/s.  The 
smaller  particles,  such  as  7.5  pm,  freeze  only  about  eight  centimeters  from  the  source,  while  the  large 
particles,  such  as  250  pm,  take  much  longer  distances  to  freeze. 

Figure  7  shows  the  evaporation  rate  of  the  particles  as  a  function  of  distance  from  the  source.  The  initial 
conditions  of  these  curves  are  the  same  as  for  Figure  6.  The  particles  with  smaller  radius  tend  to 
evaporate  faster. 
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Figure  9.  The  Effect  of  Initial  Particle  Size  on  Vcpor  Velocif 


The  Effect  of  'nitiol  Porticle  Size  C"'  Vcoor  Nun^ber  Density 


I'able  3.  Parameters  Used  for  Figures  6-12 


Initial  ■I’emperalure  (”K) 

Initial  Velocity  (m/s) 

Initial  Source  Radius  (m) 

Initial  Particle  Radius  (pm) 

Initial  Vapor  Pressure  (N/m^) 

Initial  Vapor  Density  (kg/m^) 

Initial  Volume  Fraction 

Evaporation  Coefficient 

Drag  Coefficient 

Particle  Nusselt  Number 

Vapor  Thermal  Conductivity  (W/m  ®K) 

Vapor  Specific  Heat  (J/kg  ®K) 

Heat  of  Evaporation  (J/kg) 

Heat  of  Fusion  (J/kg) 

Particle  Density  (kg/m^ 

Gas  Constant  for  Water  (J/kg  ®K) 
Stefan-Boltzman  Constant  (W/m^  “K"*) 
Emissivity  of  Particle 


343.15 

15 

7x  lO"* 

7.5,  25,  or  250 

3.119  X  lO'* 

0.19833 

0.74 

0.05 

0.9 

2.73 

0.024403 
1863.172 
2.33.382  X  10* 
3.34077  X  10* 
977.7 
461.94 
5.669  X  10-* 


0.9 


Figure  8  shows  how  particle  velocity  changes  as  a  function  of  distance  from  the  source.  For  the  large 
particles,  the  velocity  changes  very  little  from  its  initial  velocity  of  15  m/s.  However,  for  smaller  particles, 
the  velocity  changes  very  rapidly  near  the  source  and  becomes  almost  constant  far  from  the  source.  This 
phenomenon  can  be  explained  by  the  forces  that  inQuence  particle  speeds.  As  explained  earlier,  the 
vapor  velocity  increases  very  rapidly  after  ejection  from  the  source.  This  can  be  seen  in  Figure  9  where, 
in  less  than  one  half  millimeter  from  the  source,  the  vapor  velocity  reaches  sonic  speed.  ITierefore,  this 
sharp  increase  in  vapor  velocity  can  affect  the  velocity  of  the  particles,  especially  the  smaller  ones.  'Ihis 
situation  can  be  seen  in  Figure  8  where  the  7.5  (im  particle  velocity  increases  from  15  m/s  to  as  much 
as  60  m/s.  In  contrast  the  250  p.m  particle  velocity  increases  very  slightly  (  from  15  m/s  to  about  16  m/s). 
Glenn  (1969)  explained  another  factor,  which  is  worth  mentioning  here.  In  Figure  7,  it  is  clear  that  the 
evaporation  rate  of  the  smaller  particles  near  the  source  is  much  faster  than  for  the  larger  ones. 
I'hereiure,  it  is  possible  that  the  local  back  pressure  which  is  created  by  evap>oration,  and  as  earlier 
discussed,  vapor  drag,  could  inaease  the  velocity  of  the  smaller  particles. 

Figure  10  shows  that  vapor  density  decreases  very  rapidly  from  about  0.2  kg/m^  to  O.OOOOl  kg'm^  in  only 
0.25  meter  distance  from  the  source.  Figure  11  is  similar  to  Figure  10  except  that  the  vapor  density  is 
presented  in  terms  of  vapor  number  density.  Figure  12  shows  particle  number  density  variations  for  7.5, 
25,  and  250  p-m  as  a  function  of  radial  distance  r.  Figure  12  also  shows  a  representation  of  the 
experimental  measurements  of  Curry  et  al.  (1985bl  for  a  0.13  pm  diameter  particle.  As  can  be  seen  from 
Figure  12,  the  particle  number  density  calculations  of  the  present  work  are  in  good  quantitative 
agreement  with  the  experimental  work  of  Curry  el  al.  11985bJ. 


5.  SUMMARY  AND  CONCLUSIONS 


An  analytical  investigation  was  made  of  a  one-dimensional  radial  source  flow  of  water  into  a  vacuum. 
’Hie  particles  were  assumed  to  be  spherical,  with  uniform  sire  and  temperature.  The  beat  transfer 
between  particles  and  vapor  was  assumed  to  be  by  convection  and  radiation  only.  The  equations  that 
describe  this  model  are  a  sixth-order,  nonlinear,  two-point  Ijoundaiy  value  problem.  The  problem  was 
solved  by  the  "shooting/splitting"  technique.  The  results  of  this  investigation  were  compared  with  the 
other  works  in  the  literature,  and  were  generally  in  good  agreement. 

Data  similar  to  those  of  recent  shuttle  orbiter  water  dumps  into  space  were  used  to  give  some  indication 
of  water  flow  characteristics  into  a  vacuum  environment.  Since  the  results  of  the  present  work  are 
limited  to  a  one-dimensional  model,  for  a  better  understanding  of  the  problem,  the  following  goals  are 
suggested  for  future  study.  A  two-dimensional  continuum  model  of  analysis  for  a  source  flow  of  liquid 
water  jet  into  a  vacuum  (for  distances  close  to  the  source)  is  the  next  logical  step.  In  order  to  investigate 
the  complete  water  flow  structure  us  it  is  released  into  space,  kinetic  models  would  be  necessary. 
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APPENDIX  A.  DERIVATION  OF  GOVERNING  EQUATIONS 


Under  the  assumptions  stated  in  the  text,  and  using  the  control  volume  analysis,  the  system  of  equations 
governing  the  flow  are  derived  as  follows: 


A.1  CON  I  INUITY  or  VAPOR  AND  PAR'l  ICLE 

Ihe  net  flux  of  vapor  out  of  a  spherical  control  volume  is  equal  to  the  rate  at  which  mass  is  generated 
by  eva|)oration  from  the  particles.  With  this  analogy,  the  continuity  equation  of  vapor  can  be  stated  in 
integral  form  as 

JJpv  (l-“)  “  J//'”  “ 


Where  p^,  a,  u^,  and  m  are  the  vapor  density,  the  particle  volume  fraction,  vapor  velocit>’,  and  mass  rate 
evaporated  per  unit  volume  of  particles,  respectively,  n^  represents  the  unit  vector. 

Applying  Gauss'  theorem  to  the  left  side  of  the  Equation  (Al)  yields: 

JJPv  (1-a)  V  ^  Jff  (p^  (1-a)  u^I  dV  (A2) 


therefore.  Equation  (Al)  can  be  stated  as 

|J|— (Pv  (l-«)  (A3) 


In  the  spherical  awrdinate  system,  with  motion  restricted  to  the  radial  coordinate, 

-7  (Pv  (i-“)  «v)  ♦  -  IPv  (I-“)  «vl  =  ^  “  (A-i) 

.  dr  r 


or  simply: 

—  ('■^  Pv  {l'»)  u^l  =  m  U  r^  (A5) 

dr 


A-1 


By  ;in;ilogy  with  Equation  (Al).  the  continuity  equations  for  particles  can  be  stated  as  follows: 

JJpp  a  Up^  =  '|J|rn  a  dV  (A6) 


where  Pp  and  Up  are  particle  density  and  velocity,  respectively. 


Proceeding  as  before, 


JJpp  a  Up^  =  JJJA  |pp  a  «pj  dV 


(Al) 


Thus,  equation  (A6)  can  be  written  as 

C  C  C  d 


JlJ— if’p  “  “pji 


(AK) 


or,  in  one-dimensional  spherical  coordinates 

1-2 


dr 


I*"  Pp  «  “pi  =  -f”  “  »■ 


(A9) 


A.2  MOMENTUM  EOUAI  ION  FOR  VAPOR  AND  PARTICLE 

Ihe  net  force  acting  on  the  v.Tpor  within  the  control  volume,  plus  the  rate  at  which  momentum  is  added 
to  the  vapor  by  the  evaporating  particles,  must  equal  the  net  rate  of  momentum  outflow.  The  momentum 
equation  for  the  vapor  can  be  stated  as 

-f  (l-a)  dA  *  f  rf(a  f*P±  (1-a)  dV  * 

^  dx^  (AlO) 

jjjfn  a  Up.  dV  =  JJpy  (l-a)  u^.  u,.  dA 


when  fp  is  the  total  drag  force  exerted  on  the  vapor  per  unit  volume.  P  represents  the  local  pressure. 
Again,  applying  Gauss’  theorem  twice  in  Equation  (AlO),  the  result  becomes: 


A-2 


-JJJ 


ap 


(l-a)  —  -  /„  a  -  m  a 


p  -  ■  ”  *‘j>] 


dV 


(All) 


JJJ 


Pv  (l-“)  *^1  —  V  +  — [Pv  (l-“)  vl  f 

ajj  '  *  > 


Using  the  vapor  Equation  (A3)  and  substituting  it  into  Equation  (All),  the  result  then  becomes 


-JJJ 


dP 


(1-a)  —  -  /  a  -  m  a  (Up,  -  t<^.)  I  dV 
dx^  ^  n  J 

(l-«)  «Vj  —  dV 


(A12) 


Restricting  motion  along  the  radial  coordinate.  (A12)  reduces  to: 

Pv  (l-“)  \  ^  (^■“)  “  /p  =  a  (“p  -  S') 


llie  momentum  equation  for  the  particle  can  be  formulated  in  a  manner  similar  to  the  analogy  used  for 
the  vapor  momentum  equation  and  one  obtains: 


-jj  Pa  n,  dA  -JIJ 


da 


-k-p- 


dZi 


dV  “JJj  rn  a  Up  dV 


J  ^ 


(A14) 


JJ  Pp  “  V,  “pj  -i  '''> 


Using  Gauss’  theorem  and  employing  the  particle  continuity  equation  (A8),  the  momentum  equation  for 
the  particle  in  one-dimensional  radial  coordinates  becomes 

fp  “  ’V  “p  *  “  ^  “  -^P  “  ° 


A-3 


A.3  ENERGY  EQUATION  FOR  VAPOR  AND  PARTICULATE  MATreR 

From  the  first  law  of  thermodynamics,  the  rate  of  heat  addition  to  the  vapor,  plus  the  rate  at  which  work 
is  done  on  the  vapor,  plus  the  rate  of  energy  addition  due  to  the  evaporating  particulate  matter,  must 
equal  the  rate  of  energy  outflow  of  the  vapor, 

Jj  j  ?p  «  +  JJJa  /p  tip^  dV  -JJ  P  {\  -  a)  dA  ^ 

III  ^  \  \  V  =  (A16) 

/I  Pv  (1  -  a)(«v^.  ~  \  \ 


where  qp  is  the  total  rate  of  heat  transfer  per  unit  volume  by  convection  and  radiation  from  the  particles 
to  vajior.  Ay  is  enthalpy  of  the  vapor  at  particle  temperature  and  ^  .  is  internal  energy  of  vapor 

I'p  ly 

at  vapor  temperature. 


Equation  (A16)  can  be  written  as 


tp  •  /p  “p,  *  ”  \  \  \) 


dV 


(A17) 


//  Pv  (1  -  «)(V,.  ^  \  \  «v.)  ^ 


dA 


where  h  =  e  +  P/p.  Again,  using  Gauss’  theorem,  the  right  hand  side  of  equation  (A17)  becomes 


JJ  Pv  (1  -  a)(\,^  *  I  \  V  \ 


aij 


Pv  (1  -  “)  iK  *  -  \  \) 


dV 


I/I 


(A18) 


fV  (1  -  a)  u,  (V  +  1  V) 


H,.  V)  —  fpv  (1  -  “) 

‘v  i  *  *  3:^  I  J 


dV 


Using  (A3) 


Pv  (1  -  “)  \ 


m  a 


(A3) 


and  substituting  the  above  equation  into  (A18).  we  have 


III“  7  “P.  “Pi 


II/Pv  (1  ■  \  ^  " 

ai  *v 


-  V  «v)  dV 
2  ■  ' 


(A19) 


If  tlie  vapor  is  now  assumed  to  behave  as  a  perfect  gus,  so  that 

K-jc^dT 


where  is  specific  heat  of  the  vapor  and  considered  to  be  constant,  then  Equation  (A19)  can  be 
reduced  to  the  following  equation 


A-5 


(1  -  a)  (1  -  a)  u 


2  / 

tl 

V  V 


“  (9p  Wp  Up]  +  m  a  ICp^  (T^  -  TJ  ^  1  (u^'  -  «')] 


(A20) 


Now,  multiplying  the  vapor  momentum  Equation  (A13)  by  u^  and  subtracting  the  result  from  (A20) 
yields: 


(1  -  «)  C'p  \  (1  -  a) 


P  = 


«  7p  *  “  /p  (Up  -  \)  +  m  a 


Cp.  (^P  -  n)  •  ^  (s.  - 


(A21) 


The  energy’  equation  for  particles  by  similai  analogy  with  Equation  (A16)  can  be  formulated  as  follows: 
-JjJ  Vp  a  dK  -  Jjj  a  /p  up^  -  JJ  P  a  Up.  - 

Jj7  “  \  ^  «Vj  V  “ 

(A22) 


p  ( V  ^  0  “pj  ^ 


(1  -  3)(«s  ■"  -  V 


•Vi  ^ 


dA 


where  ft  and  ^  ,  and  p  are  the  internal  energy’  of  the  liquid  and  solid  phases  at  palicle 

P  p 

temperature,  and  liquid  volume  fraction  of  the  particles,  respectively.  \  is  latent  heat  of  vaporization  at 
I  p  and  defined  as 

\  \  ^  (A23) 


llien  Equation  (A22)  becomes 


A-6 


»p  *  4  “p,  ‘  ™  <*w  •  ^  ’S'i  “p(> 


dV 


J  Pp« 


3  ( V  "  7  ^  “Pi^ 


(A24) 


(1  -  3)(/^^  ^  i  «v.  Up; 


dA 


where  and  are  the  enthalpy  of  the  liquid  and  solid  phases  of  the  particles  at  particle 

‘p 

temperature,  respectively.  Using  Gauss’  theorem,  the  right  hand  side  of  Equation  (A24)  can  be  stated 
as 


IJ  Pp“ 


1 


(1  - 


\  1 


dA 


Si! 


d  X: 


Pp  “  13  (V  "  o  “Pi  "Pi^  * 

p  ^ 


(A25) 


(1  -  *  -  «p,  \)]\ 


dV 


Substituting  (A25)  into  (A24),  and  with  some  mathematical  manipulation,  (A24)  becomes 


JJJ  - 


(V,-^  •  •  5  *Pi  V 


dV  = 


sss 


“  “f.  71.  *  5  "p,  ‘ 


(A26) 


(V,.^  -  (1  -  3)  X,  *  -  lip,  Up,)  —  Pp  a  Up^ 


dV 


A-7 


Wlieri;  K(  is  heat  of  fusion  and  defined  as 


=  V,.  ”  *S, 

P  P 


Employing  the  particle  continuity  relation  (A8), 


J7J''“  { 'p  • 

p  Up  f  m  [X^  +  (1  -  3)  xj  1  dV 

JK "  “p.  ^ 

V,.  -  (1  -  P)  X,  .  I  «p^  Up_ 

P  ^ 

dV 


(A27) 


Equation  (A27)  reduces  to 


Pp  a  "p 


Cp  t'  -  x'  *  X,  3' 


2  / 

n_  a  u  u  = 

*  P  p  p 


-  [vp  *  /p  y  *  ^  (K  *  0  -  1^)  ^f) 


(A28) 


and,  as  before,  multiplying  the  particle  momentum  equation  (A15)  by  u  and  subtracting  the  result  from 
(A28),  yields 


I 


Pp  a  Up 


-  a 


ni  a  [Xy  +  (1  -  3)  Xf 


(A29) 


Equation  (A19)  is  subjected  to  the  following  constraints: 

Tp  >  r.  ;  3=1.  3'  =  0 

r  =  r  •  r'  =  0  (a30) 

p  ’  p 

Tp  <  r..  ;  3  «  3'  =  0 


» 


where  I ,  and  1’..  are  the  initial  and  triple  point  temperatures  respectively. 

Equations  (/Vi),  (A9),  (Al.^),  (Al.'i),  (A21),  and  (A29)  constitute  a  sixth  order  system  in  the  ten  unknown 
quantities  p^,  a,  u^.  m,  Up,  P,  fp,  '1'^,  qp,  and  Tp.  'Eherefore  four  constitutive  relations  are  required. 
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A.4  CONSTITUTIVE  RELATIONS 


■J'he  rate  of  mass  evaporated  per  unit  volume  of  particles  is  assumed  to  be  defined  by  the  Hertz-Knudsen 
equation: 


4  TT  €  (P^-P)  ( 


M 


1/2 


2  ft  TV 


) 


4  3 

—  TT  a 

3 


3  £ 


a 


(Po-nc 


A/ 


1/2 


2  TT  ft  Ti 


•) 


(A31) 


where  £,  R,  M,  and  are  evaporation  coefficient,  universal  gas  constant,  molar  mass,  and  vapor 
pressure,  respectively.  ’Ihe  particle  size,  a,  can  be  related  to  particle  velocity  and  volume  fraction  by 
assuming  that  no  particles  are  either  created  or  destroyed.  The  number  density  of  particles  np  (number 
of  particles  per  unit  volume)  can  be  written: 


a 


(A32) 


and  applying  the  aforementioned  conseivalion  theorem: 

— (rS%)*0  (A33) 

dr 


substituting  (A32)  into  (A33)  and  integrating.  Equation  {A33)  then  becomes 


r^aup 


=  C 


4 

—  TTCJ 


3 


3 


(A.34) 


where  the  constant  C  can  be  defined  in  terms  of  the  constants  at  the  source 

2  2 

r  aup  %a,“p, 

4343 

—Tier  — TTCT 

3  3  * 


Ihen  (A34)  becomes 

(A35) 


where  the  subscript  (*)  refers  to  conditions  at  the  initial  source  radius,  r  =  r.. 
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Equation  (A35)  can  be  written  as 
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(A42) 


^  a 


and 


?r“ 


47rCT^«,CTsB 


(A43) 


where  and  are  the  Stefan-Boltzman  constant  and  emissivity  of  the  particle.  Equation  (A3)  is 
simplified  to  the  following  equation; 

(A44) 

cr 


Equations  (A5),  (Ad),  (AI5),  (A21),and  (A29),  together  with  Equations  (A31),  (A36),  (A38),  (A42).  a^'d 
(A44)  constitute  the  complete  system.  Therefore,  they  can  be  consolidated  into  a  sixth  order  system  in 
the  unknowns  p^,  u^,  Up,  Tj,,  T^,  and  a. 
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API’KNDIX  D.  GLOSSARY 


A  Area 

j 

A|-A^  Intermediate  parameters  defined  in  Equations  (16-21) 
Specific  heat  of  vapor 
Co  Drag  coefficient 

e  Specific  internal  energy 

fjj  Drag  force  per  unit  volume 

h  Specific  enthalpy 

k  'lliermal  conductivity  of  vapor 

m  Mass  rate  evaporated  per  unit  volume  of  particles 

M  Molecular  weight 

Op  Particle  number  density 

NUp  Particle  Nusselt  number 

P  Local  pressure 

Pp  Vapor  pressure 

Convection  heat  transfer  rate  per  unit  volume  of  particle 
qp  Convection  and  radiation  heat  transfer  (qp  =  q,  +  q^) 
c|^  Radiation  heat  transfer  rate  per  unit  volume  of  particle 

R  Universal  gas  constant,  Radial  distance  in  Equation  (23) 

r  radial  coordinate 

T  Temperature 

u  Velocity 

V  Volume 
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a  Particle  volume  fraction;  particle  volume/total  volume 

3  Particle  liquid  fraction;  liquid  particle  volume/total  particle  volume 

y  Ratio  of  specific  heats 

e  Evaporation  coefficient 

c,  Emissivity  of  the  particle 

0|  Dimensionless  gas  constant  parameter 

©2  Dimensionless  initial  radius  ratio 

O3  Dimensionless  initial  pressure  ratio 

64  Dimensionless  evaporation  parameter 

6j  Dimensionless  drag  parameter 

64  Dimensionless  heat  transfer  parameter  (convection) 

67  Dimensionless  vapor  specific  heat 

0g  Dimensionless  parameter 

0,  Dimensionless  heat  of  vaporaiion 

0,g  Dimensionless  heat  of  fusion 

0|i  Dimensionless  heat  transfer  parameter  (radiation) 

Heat  of  vaporation 
Heat  of  fusion 
P  Density 

(T  Particle  radius 

Oj,,  Stefan-Boltzman  constant 
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Initial  condition 


4> 

**  Triple  point 

V  Vapor 

p  Particle 


L  Liquid 


